CRITICAL PERCOLATION IN THE PLANE. 
I. CONFORMAL INVARIANCE AND CARDY'S FORMULA. 
II. CONTINUUM SCALING LIMIT. 



STANISLAV SMIRNOV 

Abstract. We study scaling limits and conformal invariance of critical site percolation on trian- 
gular lattice. We show that some percolation-related quantities are harmonic conformal invariants, 
and calculate their values in the scaling limit. As a particular case we obtain conformal invari- 
ance of the crossing probabilities and Cardy's formula. Then we prove existence, uniqueness, and 
conformal invariance of the continuum scaling limit. 



1. Introduction 

In this paper we study critical {p = p c = \) site percolation on triangular lattice. First we 
introduce and evaluate harmonic conformal invariants, given by the limits of certain probabilities 
as mesh of the lattice tends to zero. As a corollary we obtain crossing probabilities (predicted by 
J. Cardy in [5]) and show their conformal invariance (conjecture attributed to M. Aizenman by 
R. Langlands, Ph. Pouliot, and Y. Saint- Aubin in [12]). Then we show the existence, uniqueness, 
and conformal invariance of the continuum scaling limit. Some other similar conformal invariants 
and a different approach to the scaling limit will be discussed in a subsequent paper. For the 
general background on percolation consult the book [8], for topics related to crossing probabilities, 
conformal invariance, and scaling limits see [12, 1, 7] and other references mentioned below. 

The key property of the percolation-related conformal invariants considered is that they depend 
harmonically on a parameter z (a point inside some domain $7), which allows to determine them 
uniquely from their boundary behavior, and forces them to be conformally invariant. 

To motivate search for such harmonic conformal invariants, let us first note that one can construct 
some assuming existence and conformal invariance of the percolation scaling limit. For example, 
given a simply connected domain 0, with two boundary points a and b one can define function 
h(y, z) = h a) b(y, z) to be the expected number of different clusters spanning from the boundary arc 
ab to the boundary arc ba and separating y from z (to set up this problem rigorously one actually 
needs to average numbers of "rightmost" and "leftmost" cluster boundaries separating y from z 
and taken with appropriate signs). Then by conjectured conformal invariance we can assume fi to 
be a horizontal strip of width 1. Symmetry arguments and translation invariance show that 

h(y,z) = h(y,w) + h(w,z) , 
h(y,z) = h(y + w,z + w) , 
h(y,z) =0, if y - z € iR , 
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it easily follows that 

h(y,z) = cRe{y - z) , (1) 

and hence function h(y, z) is harmonic in both y and z. So the reasoning above suggests a good 
candidate for a harmonic conformal invariant (though it does not hint at the value of c) . 

One can expect that in the discrete case similar invariants should be (almost) discrete harmonic 
functions, and one should be able to obtain a "discrete" proof of that. Indeed, modification of the 
methods described below shows that the scaling limit (as mesh of the lattice tends to zero) of the 
discrete version of h is harmonic and satisfies (1) with c = \/3/4. Note that when y and z are taken 
on the boundary of $7, the function h gives the expected number of different percolation clusters, 
crossing a conformal rectangle. 

There are other similar and more complicated invariants, which we will discuss in the next version 
of this paper. For now we will concentrate on a single harmonic conformal invariant, which is a 
complexification of crossing probabilities. Namely we will show that the probability of a point 
inside a conformal triangle to be separated by a percolation cluster spanning two sides from the 
third one is approximately a (discrete) harmonic function. Then boundary behavior considerations 
allow to reconstruct the function uniquely and force it to be conformally invariant. There seems 
to be no short motivation for this particular function to be harmonic, but the proof is simpler 
than for other analogous invariants. Also when the point z is taken on the boundary, it gives the 
crossing probabilities of conformal rectangles. Hence Cardy's formula [5] and conformal invariance 
of crossing probabilities (in the scaling limit) follow immediately. Moreover, one obtains enough 
information to prove existence, uniqueness, and conformal invariance of the scaling limit. 

There is a classical approach (due to Shizuo Kakutani, [9]) to Brownian Motion along the similar 
lines: exit probabilities for Brownian Motion started at z are harmonic functions in z with easily 
determinable boundary values. One can prove this by showing the same statement for the Random 
Walk (discrete Laplacian of the exit probability is trivially zero), and then passing to a limit 
(Kakutani works directly with Brownian Motion). However, our proof of discrete harmonicity is 
more complicated than the Random Walk analogue: instead of checking that Laplacian vanishes, we 
find harmonic conjugates (which turn out to be similar conformal invariants). This requires working 
with the first order derivatives only, while Laplacian involves the second order ones. Furthemore, 
instead of checking Cauchy-Riemann equations (which can be done), we check that contour integrals 
vanish: it requires even less precision, and allows for an easier passage to a limit. 

Interestingly, instead of a pair of harmonic conjugate functions, we get a "harmonic conjugate 
triple." It seems that 2tt/3 rotational symmetry enters in our paper not because of the specific lattice 
we consider, but rather manifests some symmetry laws characteristic to (continuum) percolation. 

After introducing harmonic conformal invariants, we obtain enough information to construct the 
continuum scaling limit, and show its uniqueness and conformal invariance. There are different 
objects which can represent the continuum scaling limit of percolation, see the discussion in [2, 1, 
3, 4], and the references therein. In the discrete setting one defines percolation clusters as maximal 
connected subgraphs of some fixed color, and the question is what sort of object will nicely represent 
the scaling limit. 

If one wants to get the full information about the "percolation configuration" in the scaling limit, 
there seem to be two approaches (with equivalent results). 

The first one is a variation of the straightforward way, when one represents a percolation con- 
figuration as a collection of compact connected subsets of the plane, representing different clusters. 
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This has to be modified, since it happens in the scaling limit that two parts of a percolation cluster 
touch without connecting (in this particular place), since they are separated by a curve of opposite 
color. However, the modification of this approach, suggested by Michael Aizenman in [1], works: 
one should represent percolation clusters by a collection of all curves contained inside clusters of 
some fixed color. Then percolation configuration is a collection of all curves inside all clusters of 
the fixed color, such collections will almost surely satisfy certain "compatibility" relations. One 
can also add for convenience the collection of all curves of the opposite color. 

Another approach is to represent a percolation configuration as a collection of "nested" closed 
curves - external perimeters of clusters (of both colors). In the discrete case such curves will be 
simple, in the scaling limit they cease to be simple but remain "non-self-traversing." These curves 
are the unique curves which correspond to crossings by both colors, and they have a canonical 
orientation (depending on the color of the "outer" side). Then the percolation culster is a space 
inside some curve 7 which is outside of all curves of opposite orientation lying inside 7. 

The philosophy of constructing the scaling limit is to deduce first existence of subsequential 
limits by compactness arguments, and to show that quantities we know determine the law of the 
scaling limit uniquely, so it does not depend on chosen subsequence. We know values of crossing 
probabilies in the scaling limit, and the crossing events generate any reasonable cr-algebra for 
percolation configurations. However it is not immediate that there is a unique law with given 
probabilities of crossing events (e.g. horizontal and vertical strips generate the Borel cr-algebra in 
the plane, but knowing measures of those does not determine a Borel measure uniquely). 

Thus it is easier to build the scaling limit gradually, so we will start with some partial aspects 
of the "full percolation configuration," which are also of independent interest, and sometimes are 
even more convinient to work with. One can think about these either as about scaling limit of 
similar objects for discrete percolation, or as about restriction of a full percolation configuration to 
a coarser cr-algebra. 

Among these objects are the outer boundary of a percolation cluster, the external perimeter 
of a percolation cluster, and a lamination of a domain (telling which points on the boundary are 
connected by clusters inside the domain). 

Setup and notation. We study critical site percolation on the triangular lattice with mesh 5. For 
general background, consult G. Grimmett's and H. Kesten's books [8, 11]. Vertices are colored in 
two colors, say blue and yellow, independently with equal probability p = p c = \- By H. Kesten's 
[10] this probability is critical for the site percolation on the triangular lattice, but it is interesting 
to note that we do not use this fact. By a blue simple path (in a particular percolation configuration) 
we mean a sequence {bj} of vertices colored blue, with adjacent bj and bj + ±. We say that a path 
goes or connects to some set, if the dual lattice hexagon around the last (or the first, depending 
on the context) vertex intersects this set. We will often identify this sequence of vertices with the 
corresponding curve - i.e. the broken line &1&2 

For brevity we denote by r the cube root of unity: r := exp (2m/ 3), and assume that one of the 
lattice directions is parallel to the real axis. 

For points a, b (or prime ends) on the boundary of a simply connected domain Q we will call 
"arc ab, " the part of <9f2 between the points a and b (in counterclockwise direction and including 
the endpoints). It need not be an arc in the strict sense, though. Note that our considerations do 
not differ whether one considers points or prime ends. By [ab] we denote the interval joining a and 
b. 
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We denote by H the space of all Holder curves up to parameterization, endowed with the uniform 
metric (innmum of sup||7i(t) — 72(i)||oo over all possible Holder parameterizations of two curves 
by the interval [0, 1]), and by B the corresponding Borel cr-algebra. 

Note that for any M the set Hm of curves admitting Holder parameterization with norm at most 
M is compact in 7i (in the uniform topology). 

Acknowledgments. I am grateful to Lennart Carleson for getting me interested in the subject of 
critical percolation, and for his encouragement and advice during this project. For some time I 
attempted to follow the approach suggested by him and Peter Jones (to establish Cardy's formula for 
equilateral triangles first, using local geometry of triangular lattice, global geometry of an equilateral 
triangle, and a study of critical cases), so some similar elements can be found in the present paper. I 
would also like to thank M. Benedicks, I. Binder, L. Carleson, P. Jones, K. Johansson, N. Makarov, 
O. Schramm, M. Sodin, and W. Werner for valuable comments. 



Consider a simply connected domain Q and three accessible boundary points (or prime ends), 
labeled counterclockwise a(l), a(r), a(r 2 ). 

If the domain has a smooth boundary, there are harmonic functions 



where u is the counterclockwise-pointing unit tangent to dQ. Moreover, functions h±,h T ,h T 2 in 
some sense form a "harmonic conjugate triple," see the discussion below. 

As Wendelin Werner pointed out to us, function h a is also given by the probability that Brownian 
Motion started at z and reflected on the arcs a(a)a(ra), a(T 2 a)a(a) at /rac7r3-angle pointing 
towards a(a) hits a(a>) before the arc a(Ta)a(T 2 a). 

Problem (2) is conformally invariant, so we can alternatively pose it (and this works even if f2 has 
non-smooth boundary) by conformally transferring it from f2 to some nice domain. The solution 
of the problem (2) takes on a particularly nice form for an equilateral triangle Q' with vertices 
a' (1) , a' (t) , a' (t 2 ) . Then the corresponding functions become linear functions h' a , which are equal 
to 1 at one of the vertices and zero at the opposite side, i.e. are properly rescaled distances to the 
sides of the triangle. For a general simply connected domain Q, there is a unique conformal map ip 
of Q to the triangle Q', sending a(a) to a'{a). Therefore we remark that 

Remark 1. h a can be defined by h a := h' a o ip. 



2. Harmonic conformal invariants 



h (a(a),a(Ta),a(r 2 a), z) = h a (z) = h a ^(z) , ae{l,r,r 2 } , 



which are the unique solutions of the following mixed Dirichlet-Neumann problem: 

h a = 1 at a(a) , h a = on the arc a(ra)a(T 2 a) 
d 

— — -h a = on the arc a(a)a(Ta) 
d(Tv) 

d 

— 2~\^ a = on the arc a(r 2 a)a(a) 




(2) 



It is also easy to see directly (or by mapping to an equilateral triangle) that for a halfplane the 
solution of (2) is a hyper geometric function. 
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For a triangular lattice with mesh S, we define an event Q a (z), a G {l,r, r 2 } , z G f2, as an 
occurrence of a blue simple path going from the arc o(a)o(ra) to the arc a(r 2 a)a(a), and separating 
z from the arc a(ra)a(r 2 a). Then we define functions 

H n ' 5 (a(a),a(Ta),a(T 2 a),z) = H s a (z) = H a (z) , a€{l,r,r 2 } , zGQ, 

to be the probabilities of Q a (z)'s. Note that i? Q 's are constant in each of the lattice triangles, 
so we will mainly consider their values at the vertices of the dual hexagonal lattice (= centers of 
triangles). We emphasize again, that we consider only simple paths, neglecting therefore the points 
separated by the "dangling ends." Otherwise the functions would have the same boundary values 
as H^s but would be strictly bigger inside the domain (and therefore not harmonic). 

Theorem 1. As 5 — ► 0, functions converge uniformly in f2 to functions h a . 

The statement above is a bit unrigorous, since the functions are defined only on discrete 
lattice depending on 5. 

Since the problem (2) is conformally invariant, we conclude that 

Corollary 1. The limit of is a conformal invariant of the points a(l), a(r), a(r 2 ), z and the 
domain Q. 

Consider a conformal rectangle, i.e. a simply connected domain Q with four points (or prime 
ends) a(l), a(r), a{r 2 ), x on the boundary, labeled counterclockwise. Then the crossing probability, 
i.e. the probability of having a blue cluster connecting the arc xa(l) to the arc a{r)a{T 2 ) is equal 
to H^ 2 (x) (or rather to the boundary value of H 5 t2 at x). Therefore 

Corollary 2. As 5 — > 0, the crossing probability tends to h T 2(x), and hence is conformally invari- 
ant. 

For a particular "nice" domain, like a rectangle or a half plane, it is easy to find h a 's and then 
check that crossing probabilities satisfy the Cardy's formula [5]. This is particularly easy for an 
equilateral triangle, where hps are linear: 

Corollary 3 (Cardy's formula in Carleson's form). Consider an equilateral triangle abc with side 
length one and a point x G [be]. As 5 — > 0, the probability of a blue crossing from [xb] to [ac], tends 
to \xb\. 

Sometimes it is more convenient to work with a different formulation of the results above. For 
a discrete percolation configuration in conformal triangle abc there is always a unique vertex w on 
the arc be with neighbors accessible by a blue crossing from ac and a yellow crossing from ab (it is 
just the endpoint of a blue crossing from ac to be, closest to ab, or equivalently of a yellow crossing 
from ab to be, closest to ac). It follows immediately from the discussion above that 

Corollary 4 (Cardy-Carleson law). The law of w converges as S — > to a distribution, which is a 
conformal invariant of configuration (Q,a,b,c). For an equilateral triangle abc this distribution is 
uniform on the side be. 

Remark 2. Note also, that the probability of a point z G fi to be contained between the two crossings 
tends to 1 — h(,(z) — h c (z) = h a (z). 

One can also ask about the speed of convergence in the statements above. 
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Remark 3. One can check that the speed of convergence in the Theorem 1 and its Corollaries above 
is 0(5 £ ) for some e > 0, and then obtain values of crossing probabilities for the discrete case up to 
such error. After establishing values of exponents for percolation one can get e = 2/3. 

The remaining part of this Section is devoted to the proof of Theorem 1. 

2.1. Proof of Theorem 1. Take e {l,r,r 2 }. If z is a center of some lattice triangle, and 
z + r\ is the center of one of the adjacent triangles, denote by Pp(z,r]) probability of the event 
Qp(z + Qp{ z )- Then the discrete derivative of Hp can be written as the following difference of 
the probabilities: 

^Hfi(z) ■= Hp(z + n)-H p {z) = P p {z, rj) - Pp(z + r,, -r,) . (3) 

Remark 4. Since the discrete derivative is expected to have order x S, and is a difference of the 
two probabilities on the right hand side, one might be tempted to assume that the latter also have 
order x 8. Surprisingly, they have a much bigger order x 5 2 ^ 3 , so there is a non-trivial cancelation. 

Lemma 2.1 (2-7r/3-Cauchy-Riemann equations). Let z be a center of some triangle, and n be a 
vector from z to the center of one of the adjacent triangles. Then for any f3 G {l, r, r 2 } 

Pp(z, V ) = P T p(z,Tn) . (4) 

Remark 5. Lemma 2.1 holds not only for triangular lattice, but for any graph which is a triangu- 
lation. 

Remark 6. Evaluating one more discrete derivative (by shifting the boundary of the domain, rather 
than the point), one can check that for nearby points z' and z" 

P {z',n>) = P {z",n>) + 0(«5 1+£ ) . (5) 

Together with the Lemma 2.1 this gives 

^ H p{z) = Pp(z, rf) - P p (z + r), -ri) = P T p(z, rrj) - P T p(z + rj, -rn) 

( = } P rP {z, t V ) - P Tfi {z + tt], -tt,) + 0(5 1+s ) = Q^H T p(z) + 0(5 1+£ ) . 

Therefore the functions H\,H T ,H T 2 satisfy discrete Cauchy-Riemann equations (or rather their 
version for triples) up to 5 s (actually even up to 5 2 ^). But in passing to a scaling limit it is better 
to use global manifestations of analyticity, so we work with vanishing contour integrals instead. 

Proof: Name the vertices of the lattice triangle which contains z by the letters X, Y, Z starting 
with the one opposite to z + n and going counterclockwise. For the event Q' := Qp(z + v) \ Qp{ z ) 
to occur, the closest to the arc a(T/3)a(T 2 [3) blue simple path 7 going from the arc a{(3)a{rf3) to the 
arc a(r 2 /3)o(/3) should separate z from z + rj. Firstly this means that there are two disjoint blue 
paths ("halves" of 7), which go from Y and Z to the arcs a(r 2 /3)a(/3) and a(/3)a(r/3) respectively. 
Secondly, the vertex X is colored yellow (otherwise we can include it into the path 7), and is joined 
by a simple yellow path to a(T(3)a(T 2 [3). 

So we conclude that the event Q' can be described as an occurrence of three disjoint simple paths, 
joining X, Y , and Z to the arcs a(r/3)a(r 2 /3), a(r 2 /3)a(/3), and a(/3)a(r/3), and colored yellow, blue, 
and blue correspondingly. 
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But in the latter description we can easily change the colors of the paths while preserving the 
probability of Q'. In fact, for a given configuration we can chose the "counterclockwise- most" 
yellow path from X to the arc a(r f3)a{r 2 (3) , and the "clockwise-most" blue path from Y to the arc 
a{r 2 (3)a{(3). Denote by Q' the union of these two paths and the part of f2 between them, containing 
a{r 2 (3). The existence of a blue path from Z to a((3)a(r(3) depends only on the coloring of O \ fl'. 
Thus if we condition by the coloring of f2', the probability of the existence of a blue path from Z to 
a(/3)a(r/3) is the same as the probability of such a yellow path, since we can invert colors in fi \ £l ; 
and use that p = \. 

Taking expectation over all possible configurations of Q', we deduce, that the event Q' has the 
same probability as the occurrence of three disjoint simple paths, joining X, Y, and Z to the arcs 
a(r P)a(r 2 (3) , a(/3)a(r 2 /3), and a(/3)a(r/?), and colored yellow, blue, and yellow correspondingly. 
But if one inverts colors in all of Q (which preserves probabilities), this is the description of the 
event Q T p(z + rrj) \ Q T p(z), and we proved the Lemma. □ 

Lemma 2.2 (Holder norm estimates). There are constants e and C depending on the domain 
only, such that Hp has e-Holder norm at most C. The boundary values of Hp are zero on the arc 
a(r f3)a{r 2 (3) and tend to 1 at the point a((3) as 5 — > 0. 

Remark 7. It is also easy to show that on dfl 

H x + H T + H T 2 -► 1 . 

Proof: Clearly, Hp takes values in [0, 1]. To prove that it is Holder, we need to show that 

\Hp{z)-Hp(z')\ < C'\z-z'\ £ , (6) 

and it is sufficient to do so for pairs z, z' well away from one of the boundary arcs, say a(a)a(ra). 
The difference on the left hand side of (6) is equal to 

Hp(z) - Hp(z') = P (Qp(z) \ Qp(z')) - P (Qp(z') \ Qp(z)) , 

so it is enough to estimate the terms on the right hand side of the above equation. But for either 
of the corresponding events to occur, the interval \zz'\ should be joint to the arcs a(T/3)a(T 2 (3), 
a((3)a(r 2 (3), and a(/3)a(r/3), by the yellow, blue, and blue clusters correspondingly - the reasoning 
is similar to that of the Lemma 2.1. Particularly, there is a monochrome cluster connecting [zz 1 ] 
to a(a)a(Tct). 

But [zz 1 ] can be separated from a(a)a(Ta) by | log \z — z'\\ — c disjoint discrete annuli of fixed 
moduli. By the Lemma 4.1 probability of the existence of a monochrome cluster traversing such an 
annulus is bounded from above by some q < 1, regardless of size. So the occurrence of a monochrome 
cluster connecting [zz'] to a(l)a(r), implies a simultaneous occurrence of | log \z—z'\\—c independent 
events of probability at most q < 1, and hence we infer 

\Hp(z) - Hp(z')\ < 2 q\ lo ^- z 'W- c = 2q- c \z - z'\ ll ° eql , 

and the desired estimate (6) follows. 

Note that as a particular consequence of this estimate for the case z' = z + rj we obtain 

Pp(z, V ) < C'5 S . (7) 

To prove the Lemma it remains to check the boundary values of Hp. Points on the arc 
a(r (3)a(r 2 (3) cannot be separated from it by a blue path, so boundary values of Hp on this arc 
are identically zero. The boundary value of Hp at the point a((3) tends to 1 as 5 — > also because 
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of the Lemma 4.1. Indeed, there are | log <5| — c disjoint discrete annuli of fixed shape around a(f3), 
each of them contains a blue cluster circumventing it independently with probability bounded from 
below by some p > 0. Thus probability of a{f3) being separated from the arc a(T(5)a(T 2 j3) is at 
least 1 — (1 — p)\ lo s <5 l- c ) which tends to 1 as 5 — ► 0. This concludes the Proof. □ 
Take some equilateral triangular contour V with vertices in the centers of the lattice triangles and 
with bottom side parallel to the real axis. Denote vertices of V by x(1),x(t),x(t 2 ) counterclockwise, 
starting with the top one. For a function H(z) we define the discrete contour integral by 

j> H{z)dz := 5 H ( z ) +5t Yl H ^ + 6t2 Yl H ^ ' 

^ z£x(t)x(t 2 ) zex(r 2 )x(l) zex(l)x(r) 

where sums are taken over centers z of lattice triangles, lying in the corresponding intervals. 

Lemma 2.3 (Contour integrals vanish). For any equilateral triangular contour T C O of length i 
with vertices in the centers of the lattice triangles, with bottom side parallel to the real axis, and 
any (3 6 {l,r, r 2 } one has 



j> H 5 p {z)dz = j> ^H 5 Tf3 (z)dz + 0(1 5 £ ) . 



Proof: Color all triangles in chess-board fashion, so that triangles with centers on T are colored 
black. Denote by B the set of all centers of black triangles, lying on or inside T, and by W the set 
of all centers of white triangles, lying inside V. 

Fix some a G {l,r,r 2 }. Take n to be of length 5/^/3 collinear with e 7 ™/ 6 (x(t 2 o) — x(ra)) and 
denote rj' := e 7 ™/ 3 ??. We can write using (3) and (4): 

£ (H p (z + r,)-H p (z)) ® (Pp(z,v)-Pp(z + ri,-v)) 

z£B\x(a)x(r 2 a) z£B\x(a)x(r 2 a) 



(4) 



(*) 



(3) 



(8) 



(P T/3 (z,Tr)) - P T p(z + TI,-TTl)) 

z£B\x(a)x(r 2 a) 

(Prfi(z, TV) ~ Pr?(z + TV, -T V )) + O^' 1 ) 

z£B\x(a)x(ra) 

(H T p(z + ttj) - H Tf ](z)) + O^ 1 ) . 

zeB\x(a)x(ra) 

In the identity (*) we also used that two sides differ by at most 0(£5~ 1 ) terms (number of vertices 
on T), which are of the order 0(5 6 ) by (7). Similarly one shows that 

{Hp(z + r/) - H p {z)) = ( H M* + rrf) - H T/3 (z)) + O^' 1 ) . (9) 
zew zew 
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Now, combining (8) and (9) we can write (using "telescoping sums") 
£ H p {z)- £ H p (z) = 

z£x(a)x(r 2 a) z&x(ra)x(r 2 a) 

(Hp{z + rj)- H p {z)) + Y, ( H p( z + rr/) - Hp(z)) 

z<=B\x(a)x(T 2 a) zeW 

Y (H Tf5 {z + r V ) - H rP {z)) + Y {H T p(z + tt{) - H T p(z)) + O^ 6 " 1 ) 

z&B\x(a)x(ra) z&A) 
zex(a)x(ra) z£x(a)x(r 2 a) 

So for any a € {l,r, r 2 } one has 

z£a:(ci)i(T 2 a) z€x(ra)x(r 2 a) zex(a)x(ra) z^x(a)x(r 2 a) (10) 

Summing three copies of equation (10) above with coefficients and different values of a plugged in: 

— -(equation 10)| a= i — % ^— (equation 10)|q, =t + - (equation 10)| a=T 2 , 

we arrive at the desired identity and prove the Lemma. □ 
By Lemma 2.2 the Holder norms of functions H^, a € {l,r, r 2 } are uniformly bounded, hence 
from any sequence of such functions with 5 — > one can chose a uniformly converging subsequence. 
Therefore to show that the functions converge uniformly, as S — > 0, to the functions h a and 
prove Theorem 1, it is sufficient prove the following 

8 

Lemma 2.4. Assume that for some subsequence Sj — > the functions H a 3 converge uniformly in 
n to soTue functions f a . Then f a — h a . 

8 8 

Proof: Clearly the discrete contour integrals § Ha converge to the usual contour integrals § fp. 
Then using Lemma 2.3 we conclude (one may need to shift contour by < 5 so that its vertices 
are in the centers of lattice triangles, but this perturbation does not affect the limit), that for any 
equilateral triangular contour rcfiof length £, with bottom side parallel to the real axis, and any 
(3 € {l, r, t 2 } one has 

fp(z)dz = I -f T p(z)dz . (11) 
r Jr T 

Take any a G {l,r, r 2 }. Adding two copies of (11) with coefficients and different values of [3 
plugged in: 

(equation 11) \p =a - Q + ' (equation ll)\p =Ta , 



we obtain 



/aW + ^(U)-/^aW)) <'■■ (12) 
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By Morera's theorem we deduce that the function in (12) is analytic, and hence for any a £ { 1, r, r 2 } 
the function f a is harmonic with harmonic conjugate (f Ta {z) — f T 2 a (z)). It follows that for any 
a G {l,r, r 2 } and any unit vector rj 

d^ fa = Wv) fra ' (13) 

Hence 

d d d 

— — -f a = -7—f T 2 a = — = for points on the arc a(a)a(Ta) , 
o(tv) ov ov 

Odd 

-fa = tt; r/ra = ^7 r0 = for points on the arc a(r 2 a)a(a) . 



d(-T 2 v) Ja ~ d{-u) JTa ~ d{-v) 

The identities above should be understood as holding in the limit aszgfi tends to the correspond- 
ing boundary arc. Also by the Lemma 2.2 the boundary values of f a are equal to zero on the arc 
a (tq) a (t 2 o) and to 1 at the point a (a). 

Summing it up, we conclude that / Q 's satisfy the mixed Dirichlet-Neumann problem (2), which 
has a unique solution. Thus f a = h a , and we have proven the Lemma and the Theorem. □ 



3. Continuum scaling limits 

Consider some domain £1 with three points (or prime ends) a, b, and c on the boundary, named 
counterclockwise. For a discrete percolation configuration we can define "the lowest blue crossing," 
y.c. _ ryLc. _ ^ c -( 0) fr 5 c ) 5 as the closest to the arc ab, simple blue curve going from the arc ca to the 
arc be. For simplicity, if no such crossing exists, we assume it to consist of one vertex, closest to c. 
This does not influence our reasoning since by the Lemma 4.1 the probability of the existence of a 
genuine crossing tends to 1 as 5 — > 0. Similarly, we define "the highest yellow crossing," 7y c '(c, a, b), 
as the closest to the arc ca, simple yellow curve going from the arc be to the arc ab. 

These curves end at neighboring vertices on the arc be, so we can join them to obtain a simple 
curve 7° A = jf b - = 7y'g(a, b, c), which starts on the arc ab and ends on the arc ca, touching the 
arc be in the middle. The law of jg- b - is a probability measure ^f b ', supported on the space H of 
Holder curves (we identify a discrete curve on the ^-lattice with the corresponding broken line). 

Theorem 2 (Outer boundary). As 5 —> 0, the law of 7 ' 6 ' converges to a law /j, ob ' on Holder 
self-avoiding paths from ab to ac, touching be at one point, which separates them into two simple 
"halves." The law /j, ob ' is a conformal invariant of the configuration (Q,a,b,c). 

Note that in the scaling limit two halves of the j°' b ' are the outer boundaries (inside Q) of 
all yellow percolation clusters intersecting the arc ab and viewed from c, and all blue percolation 
clusters intersecting the arc ca and viewed from b. 

In the proof we automatically obtain, that the law of r y°' b ' is CCI (Completely Conformally 
Invariant in the language of [17, 19]), and points where 7°' b ' touches the three arcs satisfy the 
Cardy-Carleson law (Remark 2 also gives that h a (c) is the probability of a point to be separated by 
y>.&. f rom b anc i c ). These two properties determine the law uniquely (cf. [19]). Hence using [13, 19] 
we conclude that the law of 7°' fc ' coincides with two other laws, enjoying the same properties: with 
the laws of the hull of Oded Schramm's chordal SLE 6 (started at a, aimed at some point in be 
and stopped upon hitting be) or hull of the reflected Brownian motion (started at a, reflected on 
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ab and ac at |-angle pointing towards be, and stopped upon hitting be), which was considered by 
Wendelin Werner in [19]. 

Corollary 5. The law fi°' b ' coincides with the law of the boundary of the hull of chordal SLEq (or 
reflected Brownian motion) started at a and stopped upon hitting be. 

By the work [16] of Greg Lawler, Oded Schramm, and Wendelin Werner we know the dimension 
of the outer boundary of Brownian motion, so we deduce 

Corollary 6. The curves r y°- b - and j Lc - (and hence the outer boundary of a percolation cluster) 
have Hausdorff dimension | almost surely. 

Remark 8. Lowest crossing depends only on the area "below" it, and all our considerations are 
stable under perturbations of domains, since Cardy's formula is. So one can continue working with 
the area "above" the lowest crossing, finding crossings for some configurations there. Using an 
inductive procedure it is not difficult to construct in this way scaling limits for laminations and 
backbones (and then in a similar inductive way pass to the full percolation configuration). 

Consider now a domain f2 with two boundary points a and b. For any percolation configuration 
there is a unique curve 7 e - p - along the edges of the dual hexagonal lattice, which goes from a to 
b separating the blue clusters intersecting the arc ab from the yellow clusters intersecting the arc 
ba. This curve is called the "exploration process" (one can actually introduce time) or "external 
perimeter" (of all blue clusters intersecting ab, as viewed from ba) . 

We say that a path is self-avoiding if it does not have "transversal self intersections" (maybe a 
more appropriate term would be "non-self-traversing"). 

Theorem 3 (External perimeter). As 5 — > 0, the law of the exploration process j e ' p ' converges 
to a law /U e ' p ' on Holder self-avoiding paths from a to b. The law is a conformal invariant of the 
configuration (il,a,b). 

It is possible to prove (using known crossing exponent for "5 arms") that the curve 7 e ' p ' is almost 
surely not simple. In proving that the subsequential limit of the laws fig' p ' is uniquely determined, 
we only use its properties (locality and Cardy's formula) valid also for Oded Schramm's SLE@ 
process by [13, 18], and so we arrive at 

Corollary 7. The law fi e ' p ' coincides with that of the Schramm's chordal SLEq started at a and 
aiming at b. 

As discussed above, the discrete percolation configuration in the whole plane can be represented 
by a collection of all external perimeters: (as a collection of nested pairwise disjoint oriented simple 
loops). It is not difficult (though a bit technical) to obtain from the results above the following 

Theorem 4 (Collection of external perimeters). As 5 — > 0, the law of the collection of all external 
perimeters converges to a law fi on collections of mutually- and self- avoiding oriented nested Holder 
loops. This law (and its restriction to any domain) is conformally invariant. 

We sketch the proof below. In a subsequent paper we intend to give a different, perhaps more 
conceptual, proof of this Theorem. 
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3.1. Proof for the outer boundary. By the Lemma 4.2 (which comes from the Theorem A.l in 
[3] by M. Aizenman and A. Burchard) the family of measures {fi>f b '} s is weakly precompact in TL. 
Indeed, liniM^oo Vs^'CH \ %) = uniformly in S, whereas the sets Hm are compact. 

Thus one can chose a sequence Sj — ► so that the subsequence j/^ 6 '} converges weakly to 

a measure [i o b \ It is sufficient to show that the latter measure is independent of the chosen 
subsequence. For the rest of the section we fix the sequence {6j}. 

Since the property of being self-avoiding is preserved in the limit, the measure is supported 
on self- avoiding Holder curves, starting on the arc ab, touching the arc be and ending on the arc 
ac. The two "halves" of the curve 7°' fl ' will almost surely be simple (However, the whole curve 
need not be). Indeed, the probability of a simple curve becoming non-simple in the scaling limit 
is zero: it would imply six "arms," which happens with probability zero by Lemma 4.3. Similarly 
one shows, that almost surely the curve 7 ' 6 ' touches the arc be at a unique point (otherwise there 
would be three "arms" going to a boundary point, which happens with probability zero). Denote 
by H ob ' the set of the Holder curves satisfying the properties above. 

For two disjoint simple curves r] (going from p G ab to q 6 be and denote pq below) and rf 
(going from r G be to s G ac and denoted rs below), let A„„> be the event of r y°' h ' lying completely 
inside the conformal pentagon apqrs. It is easy to see that the event A7,?/ is equivalent to the a 
yellow crossing from ap to pqrs within apqrs, which is closest to as, ending inside qr. Since this 
probability is given by the Cardy-Carleson law, it is uniquely determined and so is the probability 
of A v tf. 

It remains to notice that the events A ri , q i generate (by disjoint unions and complements) the 
restriction of the a-algebra B to the set H oh ' , which supports the measure fi°- b \ Hence the measure 
H° b - is independent of the subsequence chosen. 

Note that we automatically obtain that fj, oA is conformally invariant, since so is the Cardy- 
Carleson law. The Corollaries follow, since the properties we used are satisfied by the hull of 
reflected Brownian motion (see [19]) and by the hull of the Schramm's SLEg (see [13]). 

3.2. Proof for the external perimeter. By the Lemma 4.2 the family of measures {^s P }s * s 
weakly precompact in H. Thus one can choose a sequence Sj — > so that the subsequence j 

converges weakly to a measure [i e ' v ' . Again, it is sufficient to show that the latter measure is 
independent of the chosen subsequence. For the rest of the section we fix the sequence {Sj}. 

Since such properties are preserved in the limit, the measure n e - p - is supported on self-avoiding 
Holder curves, starting at a, and terminating at b. When 7 e ' p '(i) is some parametrization of such 
a curve by the interval [0, 1], denote by H{t) its hull at time t, i.e. the closure of Q \ Q(t). Here 
f2(t) denotes the component of connectivity of O \ 7 e,p [0,t], containing the point b. Then the hull 
grows, i.e. it is strictly larger for larger values of t, meaning that the endpoint of 7 e ' p '[0, t] is always 
"visible" from b. Indeed, otherwise the external perimeter would enter a zero width fjord, implying 
eight "arms," which happens with probability zero by Lemma 4.3. 

Denote by H e ' p ' the set of the Holder curves satisfying the properties above. It is easier to work 
with perimeters by parameterizing them (e.g. by diameter). One can do it in the following way. 
Namely, to each 7 e H ep ' and e > we associate a broken line 7 e = (f> e (7) by the following inductive 
procedure. The first point is 7^ := 7(to) = Then 7 J+1 is the first exit of j[tj, 1] from the ball 
-B(7 J , e). If it never exits this ball, we set 7 J+1 := b and end the procedure. 
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The law of curves j t is uniquely determined. In fact it is easy to see by induction that the laws 
of tj, H(tj), and hence 7 J are uniquely determined: in fact, the law of H(tj+±) \ H(tj) coincides 
with the law /i°- b - inside 5(7j,e) \ H(tj). We use that percolation is local, and the hull depends 
only on percolation inside it (stability of Cardy's formula under perturbations of the boundary is 
also used here). 

Thus the measure 4>~ is uniquely determined. Recall that limjv/^oo A* e ' p ' ( "H \ Hm) = and 
also note that 7 e converges to 7 uniformly in Hm- Therefore, fij 1 ^) converge weakly to fi, and fj, 
is uniquely determined. 

As before, we automatically obtain that /i e - p - is conformally invariant, since so are the Cardy- 
Carleson law and [i - b \ The Corollary follows, since the properties we used are satisfied by the 
Schramm's SLEg (see [13]). 

3.3. Sketch of the proof for the full limit. Again, [3] gives enough information to obtain a 
subsequential scaling limit: a measure fx on the mentioned space, and it is sufficient to prove that 
this measure is uniquely determined. One of a few possible ways is to proceed as follows. 

For any domain Q with two boundary points a and b to the law fi corresponds the law /i e,p ' = fi^ p ' 
of external perimeter in S7, which is independent from the choice of a subsequence determining /i. 
For any particular external perimeter curve 7 e ' p ' G H e ' p '(£l), consider the connected components 
{^'} of f2 \ 7 e p -. Usual "number of arms" considerations show that almost surely no other external 
perimeters can go from one such component to another within £1 (one has to also use even more 
classical estimate for 3 arms in a half-plane). Also the law fi e ' p ' depends only on 7 ep ' (or rather on 
its "infinitesimal neighborhood"). So for any component fi' bordering the boundary on the arc a'b' 
we can do the similar considerations, retrieving uniquely the law of the external perimeter within 
Q' from a' to b' (or in the opposite direction, depending on the location of a' and b': whether it 
is ba or ab). There is a subtle point here: one has to check that in some (related to conformal 
geometry) sense law of external perimeter is stable under perturbations of a domain. Continuing by 
induction, one retrieves (in the limit) uniquely the law of all arcs of curves (from the full collection 
of external perimeters) with endpoints on the boundary of $7. Now performing a new inductive 
procedure with this law, one can obtain the full scaling limit. 

4. Appendix: Technical estimates 

We gathered in this Section some known (see [1, 3, 8]) estimates we need. 

Lemma 4.1 (Russo-Seymour- Welsh estimates). For a conformal rectangle (or an annulus) of a 
fixed shape probability of the existence of a monochrome crossing is contained in an interval \p, q] . 
The constants p > 0, q < 1 depend on the shape, but not on the size of the rectangle or mesh of the 
lattice. 

This Lemma is needed to establish Holder continuity of harmonic conformal invariants. For 
bond percolation such a statement can be found in 11.70 of G. Grimmett's book [8]. The proof is 
based on self-duality property (forcing the crossing probabilities of symmetric shapes to be Jj), so 
it equally applies to critical site percolation on triangular lattice. 

Lemma 4.2. For discrete percolation all the non-repeating paths supported on the connected clusters 
(or external perimeters of clusters) within some compact region can be simultaneously parameter- 
ized by Holder continuous functions on [0,1], whose Holder norms are (uniformly in the mesh 8) 
stochastically bounded. 
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This Lemma is used to establish precompactness of laws describing discrete percolations. It is 
stated as Theorem A.l at the end of [3] by M. Aizenman and A. Burchard for non-repeating paths, 
and the main theorems of the paper clearly apply to external perimeters as well. 

Lemma 4.3. For two concentric balls of radii r and R probability of existence for discrete percola- 
tion of 5 crossings (not all of the same color) between their boundaries is < const(r/i?) 2 uniformly 
in the mesh 5 of the lattice. Hence (e.g. by the Russo-Seymour-Welsh theory) the similar probabil- 
ity for 6 crossings is < const (r/R) 2+e , and it immediately follows that in the scaling limit almost 
surely there is no point with 6 curves (not all of the same color) incoming. 

We use this Lemma (see [1] by M. Aizenman) to show that simple curves cannot "collapse" when 
we pass to the subsequential limit. For general number of arms, the exact exponents were predicted 
by M. Aizenman, B. Duplantier, and H. Aharony in [4] (see also related papers of B. Duplantier, 
e.g. [6] and references therein). Up to now there was no direct proof, but since exponents are 
known for SLEg, it is not difficult to write a proof on the basis of the current paper and work of 
G. Lawler, O. Schramm, and W. Werner [13, 14, 15]. 

We only need a very special case of 5 arms which can be calculated directly (along with 2 or 3 
arms in the half-plane, since they have exponents 1 or 2 and some "derivative representation." ) 
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